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Abstract: We use the Lorentzian AdS/CFT prescription to find the poles of the re- 
tarded thermal Green's functions of A/" = 4 SU{N) SYM theory in the hmit of large N 
and large 't Hooft coupling. In the process, we propose a natural definition for quasi- 
normal modes in an asymptotically AdS spacetime, with boundary conditions dictated 
by the AdS/CFT correspondence. The corresponding frequencies determine the disper- 
sion laws for the quasiparticle excitations in the dual finite-temperature gauge theory. 
Correlation functions of operators dual to massive scalar, vector and gravitational per- 
turbations in a five-dimensional AdS-Schwarzschild background are considered. We find 
asymptotic formulas for quasinormal frequencies in the massive scalar and tensor cases, 
and an exact expression for vector perturbations. In the long-distance, low-frequency 
limit we recover results of the hydrodynamic approximation to thermal Yang-Mills 
theory. 
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1. Introduction 

Studies of gauge theory /gravity duality at nonzero temperature and density provide 
interesting information about both thermal Yang-Mills theory at strong coupling and 
the physics of black holes/branes in asymptotically AdS space. A simple yet sufficiently 
rich example of such a duality in four dimensions is given by the correspondence between 
thermal Af = 4 SU{N) supersymmetric Yang-Mills (SYM) theory at large N and large 
't Hooft coupling and the near horizon limit of the gravitational background created 
by a collection of parallel nonextremal three-branes. This background is equivalent 
to the one of the AdS-Schwarzschild black hole of infinitely large radius [0. 

Superconformal Yang-Mills theory at finite temperature has only one scale (the 
temperature itself), and thus it cannot be regarded as an adequate model of realistic 
thermal gauge theories. Nevertheless, the dual AdS-Schwarzschild background has 
been explored as a good laboratory for studying non-perturbative features of finite- 
temperature field theory: entropy 0, ^ and transport coefficients |^ were computed, 
and the hydrodynamic approximation was used to provide a quantitative check of the 
gauge theory/gravity duality in the absence of supersymmetry p. 

From a technical point of view, the five- dimensional AdS-Schwarzschild background 
is a difficult territory. Wave equations in this background essentially reduce to the Heun 
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differential equation, a second order ordinary differential equation with four regular 
singularities, whose generic solutions are not known explicitly. Accordingly, computing 
even the simplest thermal two-point function via AdS / CFT remains beyond the reach 
unless some simplifying approximations are made. In the absence of Lorentz symmetry 
at finite temperature, two natural parameters arise, given by to = uj/2ts:T and q = 
g/27rr, where uj is the frequency of fluctuations and q is the magnitude of the spatial 
momentum. Wave equations can then be solved in the low- (t», q ^ 1) and high- 
(t», q -C 1) temperature limits 0. 

The finite-temperature correlators of A'' = 4 SYM operators dual to scalar, vector 
and tensor fluctuations in the AdS-Schwarzschild background were computed in ||, in 
the high-temperature (hydrodynamic) approximation using the Lorentzian AdS/CFT 
prescription proposed in p and recently justified^ and generalized in implementing 



the earlier ideas of |]rT|, [|T2[. In this approach, finding poles of the retarded thermal 
correlators is equivalent to computing quasinormal frequencies of a dual perturbation in 
the AdS-Schwarzschild background (see Section ||of the present paper). In the hydro- 
dynamic approximation, the thermal R-current and the stress-energy tensor correlators 
computed from gravity exhibit poles with dispersion relations t» = t»(q) predicted by 



field theory (more precisely, by relativistic fluid mechanics [0). For example, one of 
the AdS-Schwarzschild gravitational quasinormal frequencies reads [^] 

This dispersion relation is in perfect agreement with the one found in the low frequency, 
long wavelength regime of the dual thermal field theory, where the retarded correlator 
of the appropriate components of the stress-energy tensor has a pole at 

^(g)=t'.g-^^(c+^r/)g^ (1.2) 

where for AT = 4 SYM theory e = 3 P = ?,ti'^N'^T^ /S, r] = nN'^T^/S, ( = 0, and 
Vs = 1/ is the speed of sound. 

In the same regime, finite temperature Green's functions of operators dual to the 
minimally coupled massless scalar do not have poles. In this last case, going beyond the 
high-temperature approximation reveals a presumably infinite sequence of poles in the 
complex to plane corresponding to damped thermal quasiparticle excitations of A/" = 4 
SYM plasma [^. In the low temperature limit, the poles merge, forming branch cuts 
exhibited by zero temperature correlators. 



-'^An alternative justification for the prescription was suggested in ]lO| . 
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In this paper, we generalize the work of to include R-current and stress-energy 
tensor correlators, as well as the correlators of gauge-invariant operators dual to mas- 
sive scalar modes in the AdS-Schwarzschild background. We combine analytical and 
numerical methods to compute the discrete spectrum of quasinormal frequencies and 
obtain their asymptotic behavior as well as their dependence on the spatial momentum 
q. For a vector perturbation at vanishing q we find a simple analytic solution for the 
modes, given by Heun polynomials. In this case, the spectrum is known exactly and the 
thermal R-current correlator poles are determined explicitly. Spectra of quasinormal 
frequencies similar to the one found in the vector case are also observed for massive 
scalar and gravitational perturbations. Typically, frequencies stay bounded from zero, 
and do not show up in the hydrodynamic regime as poles of the correlators in the dual 
CFT. In some cases, however, a special stand alone frequency with Ren? = appears 
whose t», q -C 1 limit coincides with the analytic expression for the diffusion pole of the 
retarded correlators computed in |^ in the hydrodynamic approximation^. 

We organize the paper as follows. In Section |^, we propose a general definition 
for quasinormal modes of various perturbations in an asymptotically AdS background, 
with boundary conditions determined by the AdS/CFT correspondence. Quasinormal 
frequencies of massive scalar fluctuations in AdS-Schwarzschild geometry are deter- 
mined in Section ^. In Section ^ we find vector quasinormal frequencies corresponding 
to the poles of thermal R-current correlators. The poles of certain components of the 
finite-temperature stress-energy tensor correlators are determined in Section ^ Our 
conclusions are presented in Section |^ . 



2. A definition of quasinormal modes in asymptotically AdS 
space 

Quasinormal modes are classical perturbations with non- vanishing damping propagat- 
ing in a given gravitational background subject to specific boundary conditions. When 
the geometry is asymptotically fiat, the choice of boundary conditions is physically well 
motivated: no classical radiation is supposed to emerge from the (future) horizon, and 
no radiation originates at spatial infinity where an observer is waiting patiently to de- 
tect an outcome of some violent gravitational event (for recent reviews and references 
on quasinormal modes in asymptotically fiat space see [|14| ). 

In the case of asymptotically AdS space, one has less basis for intuition. Quasinor- 
mal modes in the relevant geometry have been studied in many publications ||15||-PP 



^In this paper, we do not consider correlators exhibiting the sound wave pole ||]. Their treatment 
requires more complicated analysis. 
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and various boundary conditions defining quasinormal modes in asymptotically AdS 
space were suggested in the literature. At the horizon, the condition is clearly the same 
as in the asymptotically flat case (no outgoing waves). Since AdS space effectively acts 
as a confining box'^ , a natural choice of the condition at the boundary would seem 



to be a Dirichlet one, at least for a scalar perturbation. However, as noted in ||22 
for the BTZ case, and as we shall see in Section ^ of the present paper, the Dirichlet 
condition is not adequate for certain values of the mass parameter. It is also not a 
suitable condition for vector and gravitational perturbations. Another condition, used 
for perturbations of the BTZ black hole in |2^, p8|, is the vanishing flux boundary 
condition. It gives the correct BTZ quasinormal frequencies, including those cases when 
the Dirichlet condition fails. It would be interesting to understand the meaning of the 
vanishing flux condition from the AdS/CFT point of view'^ , as well as to check it 
in higher- dimensional examples. An important observation made in ||2^ was that the 
BTZ quasinormal frequencies coincide with the poles (in the complex frequency plane) 
of the retarded correlators in the boundary CFT. An explanation of this fact as being 
one of the consequences of the Lorentzian AdS/CFT prescription was provided in 
Following the logic of , here we propose a pragmatic general definition of quasinormal 
frequencies which directly follows from the Lorentzian signature AdS / CFT correspon- 
dence: 

Quasinormal frequencies of a perturbation in an asymptotically AdS space are de- 
fined as the locations in the complex frequency plane of the poles of the retarded correla- 
tor of the operators dual to that perturbation, computed using the Minkowski AdS/ CFT 
prescription of 

We stress that the implementation of this definition involves only gravity calcula- 
tions. In the following sections, we use it to find quasinormal frequencies of massive 
scalar, vector and gravitational perturbations in a five-dimensional AdS-Schwarzschild 
background. 



3. Quasinormal frequencies of massive scalar perturbations 

According to the definition given in the previous section, computing quasinormal fre- 
quencies of a massive scalar perturbation of the near-extremal black three-brane back- 
ground is equivalent to finding the poles of the retarded Green's function of gauge 
invariant operators in thermal A/" = 4 SYM dual to that perturbation. Our approach 

•^This can be seen by writing the radial part of the Klein - Gordon equation in the Schrodinger 
form and considering the corresponding effective potential 

^This condition was originally used for quantization of a scalar field in a pure AdS space in global 
coordinates l3C[ |. 
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to solving this problem will be similar to the one used in where the massless case 
was considered. 

The asymptotically AdS five-dimensional part of the metric corresponding to the 
collection of parallel black three-branes in the near-horizon limit is given by 

ds^ = ^ {-fdt^ + rfx^) + ^dr' , (3.1) 

where /(r) = 1 — rg/r^, tq being the parameter of non-extremality related to the 
Hawking temperature T = ro/vri?^. This background is dual to the A/" = 4 SU{N) 
SYM at finite temperature T in the limit N ^ oo, Qym^ ~^ 

Using the coordinate z = 1 — rg/r^ and the Fourier decomposition 

one can write the wave equation for the minimally coupled massive scalar in the back- 
ground ( pTT|) as 



[1 + (1 - zf]^'^ ro' 0fc 0fc ^ 

'^'^ z{l - z){2 - z) z-'il - z){2 - z{l-z){2-z) z{l-zY{2-z) ' 

(3.3) 

where t» = c<j/27rT, q = |A;|/27rT, m = mR/2. Eq. ( p.3|) has four regular singularities 
at z = 0, l,2,cx3, the characteristic exponents being respectively {— it»/2, it»/2}; {1 — 
A/l + m^,l + A/l +m^}; {— 1»/2, t»/2}; {0,0}. In our coordinates, the horizon is located 
at z = 0, the boundary at 2; = 1. The mass parameter m is related to the scaling 
dimension A of the operator O in the dual CFT via ^ 

fA_, AG [1,2], 
A = { (3.4) 
[A+, AG [2, 00), 

where 



A± = 2 (^1 ± VI +m^j . (3.5) 
Eq. (|3.3|) can be written in the standard form of the Heun equation. 



y" + 



7 (5 e 
- + 7 + 



z z — 1 z — 2 



, a(3z-Q 

y + ziz-i)iz-2f ='^ ^'-^^ 



^The branch A_ of scaling dimensions does not arise in A/" = 4 SYM. However, for completeness 
we treat A as a continuous variable defined in the interval A g [f,oo), where A = 1 is the scalar 
unitarity bound in d = A. 
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by making the following transformation of the dependent variable 

(j)[z) = z-'f {z - - 2)-f y{z) . (3.7) 

The parameters of the Heun equation are constrained by the relation = a+(3+l. 

In the massive scalar case they are given by 

a = /? = -^il-t^ + 2-^, 7 = l-«tt>, 6 = 3- A, e = l-tt>. (3.8) 
The "accessory parameter" Q is 

The characteristic exponents of Eq. ( |3.6| ) are (0, 1 — 7) at z = 0, and (0, 1 — 5) at 2 = 1. 
With the parameters given by Eq. (|3.8|), the exponents become (0, ivo) and (0, A — 2), 
respectively. 

At z = 0, the local series solution corresponding to the exponent is given by 

00 

7/0(2:) =5Za„(t»,q)2;", (3.10) 

n=0 

where ao = 1, ai = Q/2'y, and the coefficients a„ with n > 2 obey the three-term 
recursion relation 

an+2 + An{xx>, q) ttn+i + 5„(tt), q) a„ = , (3.11) 

where 

A.(.,q) ^ Jn + im + e + 3in + ,)HQ 
^ 2(n + 2)(n + l + 7) ' ^ ^ 

D ^ ^ {n + a){n + (3) 

2(n + 2)(n + 1 + 7) 

The local solution at 2; = 0, Eq. ( p.lO|) , is expressed as a linear combination of the two 
local solutions at 2; = 1 as 

yo{z)=Ayi{z)+By2{z), (3.14) 

where for integer A 

y^{z) = (I + ...), (3.15) 

y2{z) = l + --- + /i(n),q)yi(2)log(l-z), (3.16) 
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where eUipses denote terms of order 1 — z and higher, and A, B are the elements of 
the monodromy matrix of the Heun equation. For noninteger A the logarithmic term 
in Eq. (|3.16|) is absent. 

To compute the retarded Green's function of the operator Oa dual to the per- 
turbation (f){z), one chooses the solution of Eq. ( |3.6| ) with the exponent at z = 
(this corresponds to the incoming wave condition for (j){z) at the horizon) and proceeds 
according to [If. The correlator is then proportional to A/B and thus finding its poles 
is equivalent to finding zeros of the coefficient B. For all A G [1, oo) except A = 2 the 
latter condition is in turn equivalent to the vanishing Dirichlet boundary condition at 
z = 1. Consequently, one may look for the poles (or quasinormal frequencies) simply 
by approximating an exact solution by a finite sum, and solving the equation 



N 



l/o(l)~ J]a„(tt>,q)=0 (3.17) 



n=0 



numerically, provided the series ( ^.1(J| ) converges at \z\ = 1. Asymptotic analysis of 
the large n behavior of the coefficients a„ shows that for any value of rx>, q the series 
(|3.10|) is absolutely convergent at = 1 for A > 2 (note that with the definition ( p.4|) 



A G [1,2) is equivalent to A G (2,3]). In most cases, this approach works very well. 
In practice, however, one usually needs to evaluate a large number of terms to achieve 
a good accuracy. Another difficulty is that for A = 2 the series is logarithmically 
divergent^ . 

For integer conformal dimensions, there exists an alternative approach based on 
rapidly converging continued fractions. One may notice that whenever the second solu- 
tion 1/2 (-z) contains logarithmic terms (i.e. whenever the coefficient h{w, q) in Eq. ( p.l6| ) 
is nonzero), the condition B = is equivalent to the requirement of analyticity of the 
solution ( p.lOj ) at z = 1. This condition of analyticity translates (see [^) into the 
requirement for the spectral parameter ti?(q) to obey the transcendental continued frac- 
tion equation 

Q Bo{w,<{) fii(tt),q) S2(t»,q) 



27 Ao(t»,q)- Ai(t»,q)- A2(t»)- 
More generally, in this case the coefficients a„ obey 

ttn+i 5„(t»,q) 5„+i(i»,q) 5„+2(t») 



(3.18) 



(3.19) 



an A„(j»,q)- A„+i(i»,q)- A„+2(t»,q)- 

Thus, scalar quasinormal frequencies are the solutions \x> = o)(q) of the eigenvalue 
equation ( |3.18| ) modulo those values of vo for which the solution 2/2 (-z) of the Heun 



^Actually, one can use Eq. ( 3.17 ) to estimate the values of frequencies even in that case, but the 
accuracy is very hmited. 
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equation is free from logarithms. We call the latter "false frequencies" since even 
though they automatically appear as roots of Eq. ( |3.18|) , they do not correspond to 
poles of the thermal gauge theory correlators. 

The real and imaginary parts of the solutions to the continued fraction equation 
( |3.18| ) are shown in Figures ^ and |^ as functions of the continuously varying parameter 
A. For integer A, the false frequencies are indicated by blank ellipses. 

The values of the false frequencies can be found analytically following the standard 
procedure first described by Frobenius (see e.g. §16.33 of |^). Here we only list the 



results for the lowest values of A. The false frequencies are the solutions to the algebraic 
equations 

A = 3 t»2 _ q2 ^ , (3.20a) 
A = 4 (t»2 - q2)2 _ g ^ (3.20b) 
A = 5 (t»2 - q2)^ + 3tt>2 + q2 = , (3.20c) 
A = 6 - 4qW + 2(3 + q^)(3t»^-2qW) + (q^- 3)2 = 0, (3.20d) 

and so on. There are no false frequencies for A = 2, and the A = 1 set is equivalent 
to that of A = 3. For q = the expressions are less cumbersome, and we list more of 
them: 



A 


= 3 


i»2 = 0, 


(3.21a) 


A 


= 4 


tt>^ = 0. 


(3.21b) 


A 


= 5 


i»2(3 + tt)^) =0, 


(3.21c) 


A 


= 6 


9 + 18 to^ + to^ = , 


(3.21d) 


A 


= 7 


tt)2(252 + 63 tt)^ + tt)^) =0, 


(3.21e) 


A 


= 8 


tt)^(3024 + 168 + tt)«) = , 


(3.21f) 


A 


= 9 


i»2(72900 + 21105 + 378 + w^^) = , 


(3.21g) 


A = 


= 10 


893025 + 1803060 vo^ + 104454 to^ + 756 w^^ + ro^^ = . 


(3.21h) 



Knowing values of the false frequencies exactly allows us to check the accuracy of the 
continued fraction method. For example, to = 0, n? = ±(—3)^/'' ~ ±0.93060485910 ± 
0.93060485910 z and tx> = ±z(-3)i/^ ^ tO.93060485910 ± 0.93060485910 z are the so- 
lutions of Eq. ( p.21c| ). Solving the continued fraction equation ( |3.18| ) numerically, all 



significant figures shown above are reproduced correctly (and higher accuracy can be 
achieved, if desired). These checks give us confidence that the values of quasinormal 
frequencies reported in and in this paper are determined correctly''. 



^Nevertheless, numerical difficulties in solving the continued fraction equation remain. For large 
|t»| and/or A our algorithm suffers from instability which we were not able to overcome. 
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Another check can be made by considering the t» = 0, q = hmit, in which an 
analytic solution to Eq. ( p.6|) is available. This case is discussed in the Appendix. 

Solid lines in Figures |^ and represent solutions to Eq. ( |3.18| ) for non-integer A. 
These are also false frequencies, since in this case the non-analytic part of the solution 
is given by Eq. ( p.l5| ) rather than by Eq. ( |3.16| ) (the logarithmic term is absent in 
Eq. ( p. 161 ) for non-integer A) and thus the continued fraction equation determines zeros 
rather than poles of the correlators. The true quasinormal frequencies for non-integer 
A can be found using Eq. (|3.17|) . 

Taking all these subtleties into account, we present our results for scalar quasi- 
normal frequencies in Figures 0, | and in Table |I| for q = 0, and in Figures |^ and ^ 
for nonzero q. Frequencies appear in symmetric pairs (±Ret»„, Imt»„), as reported in 
Table 0. (We do not show the symmetric Rett)„ < branch in most of our figures.) 
We observe the following properties of the spectrum: 

• Despite the appearance, the solid lines in Figures are not straight lines. For 
zero spatial momentum q, the n-th frequency is given by 

x»n= (n+^^]{l-t) + e{n,A), n = l,2,..., (3.22) 



where 



fi»,(A) + 0(l/n'^ 
6(n,A) = <^ " ' (3.23) 

(^0(1/A"2), A^oo. 

Formula ( p.22| ) generalizes the asymptotic expression found in |2^ for the mass- 



less (A = 4) case. The asymptotic parameters tt)*(A), ai^2 can in principle be 
computed numerically, but certainly a genuine analytic asymptotic expansion 
confirming (|3.22| ) is highly desirable^ . 

The magnitude of the imaginary part of quasinormal frequencies increases with A 
increasing, in agreement with the intuitive expectation that the late time behavior 
of thermal excitations should be dominated by the small A contributions. Note 
that none of the frequencies lies in the region w <^ 1, q ^ 1- This means that none 
of the Af = 4 SYM operators dual to massive scalar fields exhibits hydrodynamic 
behavior similar to the one described in lEI, ||6|1. 



^Recently, a progress has been made in obtaining asymptotics of the quasinormal spectrum for the 
Schwarzschild black hole in asymptotically flat space [Q, |Q. The approach used in ||3^ does not 
seem to work for the Heun equation, while the one employed in is promising. Another analytic 
attempt Q (based on Ikeda's approximation |^) is not quite adequate due to the lack of a small 
parameter. 
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• For all conformal dimensions, the dependence of quasinormal frequency on q is 
qualitatively the same (see Figures | and|^). This dependence is interpreted as a 
dispersion relation for quasi-particle excitations in a strongly interacting thermal 
gauge theory. Strictly speaking, however, our results apply only in the limit of 
infinite and infinite 't Hooft coupling. 



The spectrum (|3.22|) bears resemblance to spectra of scalar quasinormal frequen- 
cies of the two asymptotically AdS backgrounds where the modes can be found 
exactly, the (2 + l)-dimensional BTZ black hole pO|,|[l^, pl| and the 4d mass- 
less topological black hole [^]. However, in Eq. (p.22|) both the real and the 
imaginary parts of Wn depend on n, whereas for the BTZ, 4d massless topologi- 
cal black hole, as well as for the Schwarzschild black hole in asymptotically flat 
space, KeWn is either independent of n or reaches a finite limit as n — oo. The 
origin and significance of this difference and its possible role in the conjectured 
relation between the quasinormal spectrum and a black hole entropy in quantum 
gravity [36| are not clear to us. The fact that scalar quasinormal frequencies are 
approximately evenly spaced with n was noticed by Horowitz and Hubeny |T 
for the AdS-Schwarzschild black hole of large but finite radius. 



The pattern of quasinormal frequencies observed for scalar modes will also manifest 
itself for vector and gravitational perturbations considered in subsequent sections. In 
that case, however, a new qualitatively different type of frequency will appear charac- 
terizing the low frequency, long wavelength behavior of a dual thermal field theory. 



4. Poles of thermal R -current correlators in = 4 SYM 

The retarded thermal Green's functions of R-currents in A/" = 4 SYM in the large 
N, large 't Hooft coupling limit were computed in ||^ in the so-called hydrodynamic 
approximation, i.e. when the frequency and the spatial momentum are much smaller 
than the temperature T. Defining the retarded Green's function in the usual way, 

G^A^,q) = -ij d'xe-''^-^e{t){[j,{x), UO)]) , (4.1) 
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one finds 



Gxx = Gyy = g \ ) (4.2a) 

- - " 8 (^n> - q2) +■■■' ^4-2^) 

where ellipses denote higher order perturbative corrections in to, q. 

The appearance of the diffusion pole in the correlation functions ( [4.2b| ) - (|4.2d|) 
in the limit of small vo, q is predicted by hydrodynamics To the next order in 
perturbation theory, the position of the pole is given by 

n,= -iq2(l + q2ln2) + --- . (4.3) 

For generic values of frequencies and momenta, however, we expect additional poles 
to appear. We also expect the momentum dependence of the hydrodynamic pole in 
Eqs. (|OBD - (|Og|) to be modified. 



The approach we are following in computing thermal R-current correlators is de- 
scribed in detail in . The correlators are determined by using the Minkowski AdS / CFT 
prescription and the relevant part of the action. 



(4.4) 



where the components Ai satisfy the 5(i Maxwell equations in the near-extremal back- 
ground (PTT|). 

It turns out that each of the components A^, and Ay satisfies an equation identical 
to the one for the minimally coupled massless scalar. We immediately conclude that 
the pole structure of the correlators and G^y is the same as the one studied in |^ . 

For the component At we have the following third-order equation 

* (1-^)/ * {i-z)r ^-^^ 
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where f{z) = z{2 — z). The components Az, At and their derivatives are related by 

A; = -^A;, (4.6a) 

Az = ^i^A;' - ^ A . (4.6b) 
tt>q to 

Imposing the "incoming wave" boundary condition at the horizon, we find A'f- = 
^-jt»/2^^^^^ where F(z) is regular at 2; = 0. The function F obeys the equation 

/ Sjl-zy-l tw\ twi3-2z) 

V (1-^)/ W 2(l-z)/ 

t»2[4-(l-2)(2-z)2l ^ q2 
H ^ —F F = 0. (4.7) 

A generic exact solution of Eq. ( [4.7|) is beyond reach. (In Eq. ( |4.7|) was solved 
perturbatively in the high-temperature (hydrodynamic) limit t» ^ 1, q ^ 1.) Here, 
however, our goal is to determine for which vo and q the retarded i?-current correlators 
computed via AdS / CFT exhibit poles. This goal can be translated into the well posed 
boundary value problem as follows. 

The Lorentzian AdS/CFT prescription suggests that the retarded two-point func- 
tions of R-currents in momentum space are obtained by differentiating the expression 

(4.8) 



16 'V «l 

with respect to At and A^ representing the boundary values of the components At, A^. 
(In writing Eq. ( |4.8| ) we used the constraint ( |4.6a| ).) Thus all nontrivial information 
about the correlators is contained in the boundary limiting value of the solution to 
Eq. ( |4.7| ) (considered functional of At and A^). 

The characteristic exponents of Eq. ( [4.7| ) at the boundary z = 1 are (0,0). It 
follows that the two local solutions aX z = 1 are given by 

F\z) = ao + ai(l -z) + 02(1 - zf + . . . , , (4.9a) 
F'\z) = F' \og{l - z) + - z) + h2{l - zf + .. . , (4.9b) 

where a^, hi are the coefficients of the Frobenius expansion. The solution F[z) regular 
at 2; = can be expressed as a linear combination of F^{z) and F^\z), 

F{z) =AF\z) + B F^\z) . (4.10) 

Then, taking the limit z — 1 in Eq. ( [4.6 b| ) we get 

Bao = voqAz + q^At . 
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Therefore, for 2; ~ 1 the solution is represented by 

A[ = ^ (t»q^ + q^A) + (wqA, + q^A*) log (1 - ^) + 0(1 - z) . (4.11) 

It then follows from Eq. ( ^4.8| ) that the correlators are proportional to (constant 
and contact terms are ignored). Specifically, 

ot^Gi^-!^^^^, ,4.m, 



Comparing with Eqs. ([4.2b| ) - ( [4.2d| ), to the leading order in w, q we have B/A = ito— q^. 
In general, one looks for the poles of the correlators by demanding the condition B/A = 
0. In other words, we should determine for which values of to and q the solution F^^ 



in Eq. ( [4.10| ) is absent or, equivalently, for which to and q the solution F{z) is analytic 



at z = 1. This requirement provides us with the necessary boundary condition at 
z = 1. The problem thus essentially reduces to the one encountered in the scalar case 
in Section ^ By changing the dependent variable to y{z) = (2 — z)^^'^F(z), Eq. ( [4.5| ) 
can be written in the standard form of the Heun equation (|3.6| ) with parameters 

^^_w(l_M)^ /3 = 2 + a, 7 = 1- it», S = l, e = 1 - t» , (4.13a) 

Q^.^-=^-!=^. (4,13b) 

and the boundary conditions requiring analyticity of the solution y{z) at both ends of 
the interval 2; G [0, 1]. As discussed in |25], the analyticity condition at 2; = 1 is satisfied 
for t» and q obeying the continued fraction equation ( p.l8| ). Note that the problem of 
"false frequencies" encountered in Section ^ does not arise here since, the exponents at 
z = 1 being a multiple root of the indicial equation, the second solution of Eq. (^4.7|) is 
unavoidably logarithmic. 

For q = 0, the solutions of Eq. ( p.l8| ) can be found analytically. Suppose that for 
some n = n* the coefficient -B„„ (to) given by Eq. ( |3.13| ) vanishes, implying a constraint 
to = Wn^. Then, as Eq. ( |3.19| ) shows, all the coefficients a„(t», q) with n > vanish, 
and the solution to our spectral problem is a polynomial of degree n*, provided that 
the algebraic equation (with a finite continued fraction on the right hand side) 

Q _ Bojw) B^jw) B^jw) Bn,-2itX>) ^n,-l(t») 

27 Ao{w)- Ai{vo)- A2{V0)- A„,_2(tt>)- A„,_i(tt>) ' 
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has t» = vOn, among its solutions. With the parameters of the Heun equation given by 
Eqs. (|4.13a|) , (^4.13b| ), this occurs for 



0, V0n = n{l-i), n = 0,l,.... (4.15) 



(Indeed, if vo and q are given by Eq. ( [4.15|) , the coefficients Bn,, Bn,-2 as well as the 



expression A„,_2 — -B„,_i/v4„,_i all vanish. Solving Eq. ( |4.14]) "backwards", one can 
see that ( |4.15|) is in fact a solution.) In this case, the solutions to Eq. ( |3.6| ) are Heun 
polynomials, easily found using (|3.11|) . The first five of them, normalized to 1 at z = 0, 
are given by 

7/0 = 1, (4.16a) 

yi = l-^z, (4.16b) 

^ 6 + 3z 3 + 9i 2 /. N 

1/2 = 1 + ' ^^-^^^^ 

51 + 21i 231 + 297i ^ II - m ^ 

W3 = 1 z + + , 4.16d 

26 260 260 ' ^ ' 

, 68 + 262 108 + 111^ 2 262 + 1179i o 917-14412, 

V4 = 1 z H z^ z-^ z^ . (4.16e) 

^ 25 50 850 6800 ^ ' 

The complementary sequence of solutions^ consists of the spectrum 

q = 0, t»„ = -n(l + i), n = 0,l,.... (4.17) 
and an infinite set of polynomials 

2/0 = 1, (4.18a) 

2/i = l-i^z, (4.18b) 

etc. Heun polynomials can be regarded as local solutions of the Heun equation valid 
simultaneously at three singularities, the characteristic exponent at each singularity 
being zero. 

We remark that the pattern m„ ~ n{l — i) had appeared in |^ and in Section |^ 



of the present paper as the conjectured large n asymptotics of the scalar quasinormal 
frequencies. The above discussion suggests that Heun polynomials can be used as a 



^Note that if y{u) is a solution of Eq. (4.5) with the characteristic exponents vi = — //*/2 at m = 1 
and h'a ~ at u = (corresponding to the spectral parameter ^ = fj,,, = iro*), then y{u) is also a 
solution with the same analyticity property and the spectral parameter fj, = fi^. Thus the poles t»„ 
come in pairs, distributed symmetrically with respect to the Imt» axis. 
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good approximation of the scalar large n solution. The exact solution (|4.15|) also adds 
support to the claim ||25| that the number of scalar quasinormal frequencies is infinite. 

For q 7^ 0, the solutions to our boundary value problem are Heun functions. The 
corresponding values of the spectral parameter t»n(q) can be found numerically by 
solving Eq. (|4.14|) . A typical distribution of poles in the complex t» plane is shown in 
Figure 0. In addition to the infinite symmetric sequence of poles familiar from the scalar 
case, there exists a special stand-alone "hydrodynamic" pole located on the negative 
imaginary axis. For this pole, the dispersion curve to = t»(c|) calculated from Eq. (|3.18 



is shown in Figure ^ together with the analytic approximation ( [4.3|) obtained in |^ in 
the limit of small q. 

The lowest ten dispersion curves generalizing the sequence ( [4.15| ) to q 7^ are 
shown in Figures p|,|TO[ Their behavior is similar to the one observed in the case of 
scalar perturbations, except for the nontrivial "roton" minimum of Ret»ri(<l) shown in 
detail in Figure O . 



5. Poles of thermal stress-energy tensor correlators. 

In this Section, we shall find the poles of the retarded Green's function for the compo- 
nents of the stress-energy tensor, 

G,,,xp{uj,q) = -i! d^xe-"'-^9{t){[T,,{x), T^M]) ■ (5-1) 



We will focus specifically on the components of T^^ whose correlators possess diffusion 
poles in the hydrodynamic regime of the theory. The procedure for computing these 
correlation functions from gravity is very similar to the one used in Section ^ for the 
R-current correlators. (The reader is also referred to for details.) In the setup 
of P], correlation functions having the diffusion pole correspond to a gravitational 
perturbation g^^, = gj^J + h^j, of the AdS-Schwarzschild background g^l^y with htx 7^ 0, 
hxz 7^ 0. It turns out |P that the correlators are essentially known once the solution to 
Eq. ( |3.6| ) with the parameters 

^ ^ _ro(l_+|) _ ^ ^ /5 = 3 + a, 7 = 1- ito, 5 = 0, e = 1 -tt>, (5.2a) 
Q.,._2_5l(i±i_!^ (5.2b) 



- 15 - 



is found. The correlators are given by 



G'te,te(tt>,q) = , (5.3a) 

= ' (5-3b) 

Gx^,a'^(tt>,q) = , (5.3c) 

where A, B are the coefficients of the connection formula for the solutions of the 
Heun equation with parameters ( |5.2a|) , ( |5.2b|) . To the lowest order in t», q we had 



B/A = iw — q^/2 p. For generic to and q, the poles of the correlators are found by 
solving the eigenvalue equation ( ^.181 ). False frequencies are described by the equation 



I»2 - c,2 _ 0. 

By examining the set of coupled equations for the perturbations hf^, h^z (Eqs. 
(6. 13a)- (6. 13c) in p) we observe that in the limit q — >• equations decouple and, 
moreover, the only nontrivial equation left coincides with the one of the minimally 
coupled massless scalar in the background (|3.1|) . We conclude that for q = the 
spectrum is identical to the one of the A = 4 scalar case (given in Table and in 



25|). Curiously, this coincidence does not seem to be obvious when comparing the 
parameters of the Heun equation in the two cases. 

For q 7^ 0, the distribution of poles in the complex to plane is qualitatively similar 
to the one shown in Figure ^ There is a "hydrodynamic" pole whose dispersion relation 
(see Figure |12D for small q is well approximated by the analytic result vo = — 2q^/2 0. 
For other poles the dependence on q is shown in Figures |13|, |1^. 



6. Conclusions 



Computing quasinormal frequencies in asymptotically AdS space using the standard 
framework of general relativity may be interesting on its own right, but when the 
computation is motivated by the AdS / CFT correspondence what one really is interested 
in is a way to compute the Lorentzian signature correlators from gravity. In many cases 
this is technically difficult or impossible, and yet even in those cases computing the poles 
of the correlators may turn out to be relatively straightforward. In this paper we found 
the poles of the retarded correlators of the thermal Af = 4 SYM theory operators dual to 
scalar, vector and gravitational perturbations in the 5d AdS-Schwarzschild background. 

Since our knowledge of the strong coupling regime of the theory obtained from the 
sources other than the AdS/CFT is very limited, the interpretation of our results is not 
obvious. We clearly see the emergence of a hydrodynamic behavior in the theory, but 
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the role of the infinite sequence of "quasi-Matsubara" frequencies Un ~ 27rTn(l — i) is 
not clear. Equally, with temperature being the only scale in the theory, it is not clear 
whether the similarities in the spectrum for perturbations of different spins have an 
underlying algebraic explanation. 

Our results were obtained in the limit of infinite and infinite 't Hooft coupling. 
The corresponding perturbative calculation in a weakly coupled gauge theory does not, 
to the best of our knowledge, exist in the literature. It would be interesting to compare 
the two regimes explicitly, as well as to compute the correction to our results appearing 
at large but finite 't Hooft coupling in the spirit of 0. 

The singularity structure of wave equations describing scalar perturbations of a 
generic black p-brane near-horizon geometry |^ suggests that the spectrum of its quasi- 
normal excitations should be similar to the one observed in this paper. Studies of such 
a spectrum may prove to be useful in the effective description of black objects using 
the language of a dual real-time thermal field theory (or quantum mechanics 
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A. Solutions of Eq. ([3761) for = 0, q = 

Consider Eq. (O) with t» = 0, q = 0. Then Q = (2 - /\/2f = a(3, and the Heun 



equation reduces to a hypergeometric one. Changing variable to f = {l — zY, we obtain 

4v{v - l)y" + 2 [2t; + (4 - A){v - 1)] y' + (^2 - y = 0, (A.l) 

whose formal solution is given by a linear combination 

y{v)=C,y,{v) + C2y2{v), (A.2) 

where 

f A A A \ , 

y,{v)= 2Fiil--,l-j;2--;v] , (A.3) 



^^^v]. (A.4) 



4' 4' 2' 



-17- 



The hypergeometric functions in Eqs. (|A.3|) , (|A.4| ) are degenerate. They are exphcitly 
represented by the series expansion 



(a, a; 2a;v) = V <^ 2V'(A; + 1) -2,p{a+k) -log {1 - v) } (l-vf , (A.5) 




where a = 1 - A/4 for (|Al3| ) and a = A/4 for (|AD, the expansion still being valid 
when a is a negative integer or zero (in which case an appropriate limit should be 
taken). However, when a = -1/2, -3/2, . . . , i.e. when A = 2{2k + 1) = 6, 10, 14, ... , 
the correct representation for ( [A.3|) is instead given by 

yiiv) = y2iv)\ogv 

+ £ ^"t^^^^TTH- I + ^ + 1/2) - 2^ik + 1/2) - ^(2A; + 1 + n) 



n=l 



2k 



+ ^{2k + 1) - ^{n + 1) + ^(1) - (i/2_fc);]r ^''"" • (^-6) 



(^-1)!(-2 A;)„ 
„=i [(1/2 -A;). 

Another special case is A = 2 which is covered by ( [A.6|) with /c = and the last sum 
omitted. 

Having found the explicit solutions, we can now use them to illustrate the reasoning 
adopted in Section ^ We notice that for A 7^ 2(2A; + 1) = 2, 6, 10, 14, . . . the solution 
( |A.2| ) does not contain logarithmic terms at f = 0, in agreement with Eq. ( p.21| ). 
Moreover, ( [A.5| ) shows that in this case one can choose integration constants Ci, C2 
to get rid of the logarithms also at f = 1. Consequently, an analytic solution in the 
interval v G [0, 1] exists, and thus n? = 0, q = must be among the solutions of the 
continued fraction equation ( p. 181) for A 7^ 2{2k + 1). This is illustrated in Figures |l],|^. 
The solutions are "false frequencies" since the absence of logarithms at f = reflects 
the property of Eq. ( A. 3 ) rather than the requirement Ci = 0. 

On the other hand, for the "exceptional" conformal dimensions A = 2(2A; + 1) = 
2, 6, 10, 14, ... , Eq. ( |A.6| ) shows that logarithms do appear in the second solution, and 
that in this case there is no nontrivial analytic solution to Eq. {TB) for to = 0, q = 0, 
again in agreement with Figures |I| and 0. 
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Table 1: Real (with it sign) and imaginary parts of the ten lowest scalar quasinormal 
frequencies for integer conformal dimensions A G [2, 10] at zero spatial momentum q. 



n 


A= 


2 


A= 


3 


A= 


= 4 


1 


±0 640759 


-0 411465 


±1 099407 


-0 879767 


±1 559726 


-1 373338 

-I- • KJ 1 KJ KJ KJ KJ 


2 


±1.618564 


-1.393310 


±2.105949 


-1.887444 


±2.584760 


-2.381785 


3 


±2.614565 


-2.391212 


±3.107772 


-2.888629 


±3.593965 


-3.384782 


4 


±3.613045 


-3.390584 


±4.108584 


-3.889040 


±4.598600 


-4.386241 


5 


±4.612274 


-4.390299 


±5.109031 


-4.889237 


±5.601338 


-5.387081 


6 


±5.611817 


-5.390141 


±6.109309 


-5.889349 


±6.603123 


-6.387619 


7 


±6.611519 


-6.390042 


±7.109497 


-6.889419 


±7.604368 


-7.387990 


8 


±7.611311 


-7.389975 


±8.109631 


-7.889468 


±8.605279 


-8.388258 


9 


±8.611159 


-8.389927 


±9.109731 


-8.889502 


±9.605971 


-9.388459 


10 


±9.611043 


-9.389892 


±10.109808 


-9.889528 


±10.606513 


-10.388616 


n 


A= 


5 


A= 


6 


A= 


= 7 


1 


±2 028589 


-1 879263 


±2 506053 


-2 390899 


±2 990770 


-2 904529 


2 


±3.061998 


-2.880185 


±3.540615 


-3.382852 


±4.021993 


-3.888969 


3 


±4.077135 


-3.882247 


±4.559507 


-4.382044 


±5.042446 


-4.884374 


4 


±5.085553 


-4.883785 


±5.571025 


-5.382599 


±6.056103 


-5.883122 


5 


±6.090832 


-5.884881 


±6.578659 


-6.383372 


±7.065667 


-6.883004 


6 


±7.094413 


-6.885679 


±7.584035 


-7.384101 


±8.072656 


-7.883288 


7 


±8.096982 


-7.886277 


±8.587998 


-8.384732 


±9.077947 


-8.883702 


8 


±9.098903 


-8.886739 


±9.591023 


-9.385268 


±10.082069 


-9.884140 


9 


±10.100388 


-9.887102 


±10.593398 


-10.385719 


±11.085357 


-10.884560 


10 


±11.101564 


-10.887395 


±11.595306 


-11.386101 


±12.088032 


-11.884947 


n 


A=8 


A= 


9 


A= 


= 10 


1 


±3.481133 


-3.418069 


±3 975640 


-3 930432 


±4 473030 


-4.441154 


2 


±4.506708 


-4.397468 


±4.994847 


-4.907305 


±5.486181 


-5.417586 


3 


±5.526820 


-5.389010 


±6.013158 


-5.895504 


±6.501733 


-6.403317 


4 


±6.541573 


-6.385462 


±7.028005 


-6.889514 


±7.515798 


-7.395035 


5 


±7.552508 


-7.384005 


±8.039695 


-7.886431 


±8.527626 


-8.390213 


6 


±8.560810 


-8.383496 


±9.048935 


-8.884850 


±9.537392 


-9.387379 


7 


±9.567269 


-9.383432 


±10.056333 


-9.884076 


±10.545457 


-10.385710 


8 


±10.572405 


-10.383580 


±11.062345 


-10.883749 


±11.552165 


-11.384742 


9 


±11.576569 


-11.383826 


±12.067302 


-11.883672 


±12.557793 


-12.384204 


10 


±12.580001 


-12.384111 


±13.071442 


-12.883738 


±13.562562 


-13.383935 
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Figure 1: Real part of the eigenfrequencies (solutions of the continued fraction equation 
at q = versus the conformal dimension A. Black dots correspond to quasinormal 
frequencies at integer values of A, while blank ellipses are the "false frequencies". The dashed 
line indicates that the sequence presumably continues to infinity. 
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Figure 2: (Minus) imaginary part of the eigenfrequencies (solutions of the continued fraction 
equation (|3.18| )) at q = versus the conformal dimension A. Black dots correspond to 
quasinormal frequencies at integer values of A, while blank ellipses are the "false frequencies". 
The dashed line indicates that the sequence presumably continues to infinity. 
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Figure 3: Ret» of the lowest eight scalar quasinormal frequencies versus the conformal 
dimension A. Dots correspond to integer conformal dimensions. 
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Figure 4: -Imt» of the lowest eight scalar quasinormal frequencies versus the conformal 
dimension A. Dots correspond to integer conformal dimensions. 
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Figure 5: Reu) of the scalar fundamental quasinormal frequency vs q for the (integer) 
conformal dimensions A G [2, 10]. The lowest curve corresponds to A = 2. 




Figure 6: — Imto of the scalar fundamental quasinormal frequency vs q for the (integer) 
conformal dimensions A G [2, 10]. The lowest curve corresponds to A = 2. 
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Figure 7: Poles of an i?-current thermal correlator in the complex to plane for q = 1. 




Figure 8: — Imro of the i?-current thermal correlator's "hydrodynamic pole" as a function 



of q. The light curve corresponds to the analytic approximation (4.3) for small q. 
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Figure 9: The ten lowest dispersion curves (Re t» vs q) for the R-current correlators 
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Figure 10: The ten lowest dispersion curves (— Imro vs q) for the R-current correlators. 



-27- 



1 


1 








1 















9 











8 

7 
6 




\ 

\ 

\ 

\ 

\ 

\ 

\ 

-Imw^\ 

q 



0.2 0.4 0.6 0.8 1 



Figure 11: Details of the lowest dispersion curve for the R-current correlators. 
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Figure 12: The dispersion curve for the tensor diffusion pole. The dashed line corresponds 
to the analytic approximation t» = — zq^/2 valid in the hydrodynamic regime q <C 1. 
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Figure 13: The ten lowest dispersion curves (Re xx> vs q) for the stress-energy tensor corre- 
lators. 
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Figure 14: The ten lowest dispersion curves (— Imt» vs q) for the stress-energy tensor cor- 
relators. 
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